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The disformal transformation of metric ^ + r(((), X)9^</><?,/0, where 0 is a scalar 

field with the kinetic energy X = preserves the Lagrangian structure of Gleyzes-Langlois- 

Piazza-Vernizzi (GLPV) theories (which is the minimum extension of Horndeski theories). In the 
presence of matter, this transformation gives rise to a kinetic-type coupling between the scalar field 4> 
and matter. We consider the Einstein frame in which the second-order action of tensor perturbations 
on the isotropic cosmological background is of the same form as that in General Relativity and 
study the role of couplings at the levels of both background and linear perturbations. We show 
that the effective gravitational potential felt by matter perturbations in the Einstein frame can be 
conveniently expressed in terms of the sum of a General Relativistic contribution and couplings 
induced by the modification of gravity. For the theories in which the transformed action belongs 
to a class of Horndeski theories, there is no anisotropic stress between two gravitational potentials 
in the Einstein frame due to a gravitational de-mixing. We propose a concrete dark energy model 
encompassing Brans-Dicke theories as well as theories with the tensor propagation speed ct different 
from 1. We clarify the correspondence between physical quantities in the Jordan/Einstein frames and 
study the evolution of gravitational potentials and matter perturbations from the matter-dominated 
epoch to today in both analytic and numerical approaches. 


I. INTRODUCTION 

The large-distance modification of gravity has been un¬ 
der active study in connection to the dark energy prob¬ 
lem Modifications of the Einstein-Hilbert term 

i?/(167rG) in the Lagrangian of General Relativity (GR), 
where R is the Ricci scalar and G is the Newton gravi¬ 
tational constant, generally give rise to a radiative scalar 
degree of freedom cj) d, dl ■ Provided that the fifth force 
mediated by this new degree of freedom is suppressed in 
the solar system through Vainshtein Q or chameleon dl 
mechanisms, the same scalar field can potentially be the 
source for the late-time cosmic acceleration. 

Horndeski theories Q are known as the most general 
scalar-tensor theories with one scalar degree of freedom 
whose equations of motion are kept up to second-order 
in time and spatial derivatives (see also Refs. Q). Many 
of dark energy models proposed in the literature- such 
as f{R) gravity m , Brans-Dicke (BD) theory [Hill, 
kinetic braidings [l3j| . and Galileons [ij, [11- belong to 
a sub-class of Horndeski theories. In the presence of ad¬ 
ditional matter, the authors in Ref. [ll derived linear 
perturbation equations of motion on the fiat Friedmann- 
Lemaitre-Robertson-Walker (FLRW) background to con¬ 
front dark energy models in the framework of Horndeski 
theories with the observations of large-scale structures, 
weak leasing, and Gosmic Microwave Background (CMB) 
(see also Refs. [itM^ ). 

In BD theories, a scalar degree of freedom (f> is coupled 
to the Ricci scalar R of the form cj)R [ 2 I. The frame in 
which matter fields are minimally coupled to the metric is 
dubbed the Jordan frame (JF). The standard interpreta¬ 
tion of measurements is usually performed in this frame. 
In the JF of BD theories, the scalar field cj) mediates a 


fifth force with matter through its gravitational inter¬ 
action with the metric. This interaction can be clearly 
seen in the Einstein frame (EF) where the Lagrangian is 
described by the Einstein-Hilbert term plus a canonical 
scalar field [2J-[^. In the EF, matter fields feel a metric 
conformally related to the EF metric g^i, of the form 
g^i, = where D(0) is a conformal factor that 

depends on </) [1,0, [131 ■ 

For the theories in which field derivatives are cou¬ 
pled to the metric, the conformal transformation can 
be extended to a more general mapping of the metric- 
dubbed the disformal transformation |28l |. In fact, the 
structure of the Lagrangian in Horndeski theories is 
preserved under the so-called disformal transformation 
-\- T{(j))dfj,(i)d,y(j), where D and T are func¬ 
tions of cj} [10 [10. In the presence of matter fields, 
the disformal transformation helps us to understand the 
physical content of mixing between cj} and matter [3ll-l38|. 

Gleyzes, Langlois, Piazza, and Vernizzi (GLPV) [10 
proposed a generalized version of Horndeski theories by 
expressing the Horndeski Lagrangian in terms of the 
Arnowitt-Deser-Misner (ADM) decomposition of space- 
time [10 with the choice of the unitary gauge. This gen¬ 
erally generates derivatives higher than second order, but 
there is only one radiative scalar degree of freedom on the 
FLRW background [U [10 . The Lagrangian structure of 
GLPV theories is preserved under the disformal transfor¬ 
mation of the form 

gf,„ = n'^icl))gfj,,,-\-T{^,X)dfj,cl)d,,(j), (1.1) 

where P depends on cj) and its kinetic energy X = 
d^cfd>^f}/2 0 . 

In the single field system, it was shown that the invari¬ 
ance of curvature perturbations C and tensor perturba- 
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tioiiS 7 ij holds under the disformal transformation (HID 
(see also Refs. [ssl IdTl - f^ for related works). For 
appropriate choices of and F, it is possible to trans¬ 
form the action to that in the EF where the second-order 
action of tensor perturbations is of the same form as that 
in GR [13, [mI . This property is useful for the com¬ 
putation of primordial scalar and tensor power spectra 
generated during inflation [i^. Since the leading-order 
tensor power spectrum in GLPV theories is proportional 
to the Hubble parameter squared in the EF, the de¬ 
tection of primordial gravitational waves can determine 
the energy scale H during inflation. 

In GLPV theories, even though matter is minimally 
coupled to gravity in the JF, there is a mixture between 
the propagation speeds of the scalar field (j) and matter 
( 39 I 1^ . This comes from a kinetic-type mixing asso¬ 
ciated with the presence of higher-order derivatives be¬ 
yond the Horndeski domain, which is weighed by a pa¬ 
rameter Oh characterizing the deviation from Horndeski 
theories [^, [SJ, [55| . The disformal transformation cu, 
which contains higher-order derivatives, is helpful to un¬ 
derstand the origin of such a kinetic mixing affecting the 
scalar and matter sound speeds (i^ . 

If the scalar degree of freedom cj) in Horndeski theo¬ 
ries is responsible for dark energy, the tensor propagation 
speed Ct is typically close to 1 during the early cosmolog¬ 
ical epoch [^. This is not the case for GLPV theories, 
in which the deyiation from ct = 1 is allowed due to the 
absence of extra conditions Horndeski theories obey. Re¬ 
cently, dark energy models with constant Ct and yary- 
ing Ct haye been proposed in the framework of GLPV 
theories. In particular, the latter proyides an interesting 
possibility of realizing weak gravity for the perturbations 
relevant to redshift-space distortions [53, • 

In GLPV theories with ct different from I, the dis¬ 
formal transformation (ED to the EF should allow us 
to understand the structure of the matter-scalar cou¬ 
plings mentioned above. For the models proposed in 
Refs. [ 5 ^ [ 5 ^ the anisotropy parameter r] = —$/4' be¬ 
tween two gravitational potentials d' and $ deviates from 
1 in the JF, but we will see that it is possible to de-mix 
the gravitational potentials and the scalar field in such a 
way that there is no anisotropic stress between 'I' and $ 
in the EF. Moreover, we will show that the effective grav¬ 
itational coupling with matter can be well understood in 
the EF due to the separation of a GR-like contribution 
and modifications arising from an- 

In this paper we obtain relations of physical quanti¬ 
ties between the JF and the EF under the disformal 
transformation (HI) and study roles of gravitational cou¬ 
plings with matter at the levels of both background and 
linear perturbations. A similar prescription was taken 
in Ref. [38l| with the disformal transformation = 
J- T{4>)d^(j3dv4>^ but our treatment is more gen¬ 
eral in that it is not restricted to the transformation 
between Horndeski theories alone. We employ the ap¬ 
proach of effective field theory of cosmological perturba¬ 
tions (sqI - I^ . allowing to encompass GLPV theories as 


a specific case. We propose a new dark energy model in 
the framework of GLPV theories, which accommodates 
models with ^ 1 as well as models based on BD the¬ 
ories (which lead to “coupled quintessence” models 
in the EF). 

Our paper is organized as follows. In Sec.|lT]we briefly 
review GLPV theories and in Sec. nil we show how the 
GLPV action is transformed under the disformal trans¬ 
formation dm in the presence of matter. In Sec. m we 
present linear perturbation equations of motion on the 
flat FLRW background in both the JF and the EF. In 
Sec. El we consider the transformation to the EF and 
discuss the matter-scalar coupling in the EF. In Sec. IVTl 
we propose a new dark energy model belonging to GLPV 
theories and study the correspondence of physical quan¬ 
tities between the JF and the EF in detail. Section IVHI 
is devoted to conclusions. 


II. GLPV THEORIES IN THE PRESENCE OF 
MATTER 

GLPV theories [s^ are the generalizations of Horn¬ 
deski theories written in terms of ADM scalar quantities 
defined below [^. We begin with the line element 

ds^ = gfj^vdx^dx'' 

= -N^dt^ + h,j{dx^ + N^dt){dx^ + N^dt ), (2.1) 

where N is the lapse function, V® is the shift vector, 
and hij is the three-dimensional spatial metric. We 
express the three-dimensional Ricci tensor on the con¬ 
stant time hyper-surfaces Et, as The 

extrinsic curvature is defined by where 

= (—iV, 0,0,0) is a unit vector orthogonal to Ej. We 
introduce a number of geometric scalar quantities, as 

K = K%, S = , 

n = n%, u = . (2.2) 

Horndeski theories, which have one scalar degree of 
freedom (j), can be reformulated by using the above geo¬ 
metric scalars with the choice of unitary gauge 

0 = ^{t ), (2.3) 

under which (j) depends on the cosmic time t alone. On 
the flat FLRW background with the scale factor a{t), the 
Lagrangian of Horndeski theories can be expressed in the 
form [68| 

L = A2{N,t)+A3{N,t)K 

+A4{N,t){K^ - S) + B4iN,t)n 

+A5(7V, t)K 3 + B^iN, t) [u - i , (2.4) 

where K 3 = and Ai, 

Bi are functions of N and t satisfying the two conditions 

[H 

A 4 = 2 AH 4 .X - ^4 , A3 = -^XB3,x, (2.5) 
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where Bi x = dBi/dX with X = In the 

unitary gauge we have X = where a dot 

represents a derivative with respect to t. Hence the de¬ 
pendence on N and t translates to that on X and 4>. 

Violation of the conditions (|2.5I) can generally give rise 
to derivatives hi ghe r than second order, but it was shown 
in Refs. 1^ 1^ that there is only one scalar propa¬ 
gating degree of freedom on the flat FLRW background. 
GLPV theories are described by the Lagrangian (12.41) 
without having the two conditions ()2.5|) . In this paper we 
focus on GLPV theories in the presence of a matter field 
described by the Lagrangian Lm- Then, we consider 
the following action 

S = J d'^xy/^L{N, K,S,nM\t) 

j ^ gZ/m(5 (2.6) 


with 

a = + TX . (3.3) 

In the JF, let us consider the linearly perturbed line 
element on the flat FLRW background, 

ds^ = -(N^ + 2A)df + 2il}\,dtdx^ 

-|-a^(t)[(l -|- 2(^')Sij ^ij\dx^dx^ , (3.4) 

where N is the background value of the lapse; A, ip, (, E 
are the scalar metric perturbations; 7 ^ is the tensor per¬ 
turbation, and the lower index “|j” denotes the covariant 
derivative with respect to the three-dimensional metric 
hij. Gomparing Eq. (EH) with Eq. dSH), we have the 
relations 2A = — hijN'^N^ — N^, ip\i = hijN^ , and 

hij = a^(t) [(1 -|- 2(p)Sij + ‘2,E\ij + 7 ^] . (3-5) 


where g is the determinant of metric and L is given 
by Eq. (12.41) . The matter field 41^ is assumed to be a 
barotropic perfect fluid, which can be modeled by a k- 
essence Lagrangian P{Y) depending on the kinetic term 
Y = g ^^of a scalar field x [Hi, [11,16^. The term 
in Eq. (12.41) can be expressed as = 3iL(2iJ^ — 
2KEI + ~ S) up to second order in the perturbations 

[^ . where E[ is the Hubble parameter defined later in 
Eq. (14.11) . Hence the Lagrangian (12.41) of GLPV theories 
depends on N, K, S, TZ, lA, and t up to linear order in the 
perturbations. 

We assume that, in the JF, the matter field 'Ym is 
minimally coupled to the metric g^m. The matter energy- 
momentum tensor following from Lm is given by 


rpfj,U _ 
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(2.7) 


We can derive the background and linear perturbation 
equations of motion by varying the action (12.61) up to 
first and second orders in the perturbations, respectively 


Introducing the scale factor dpt) in the transformed 
frame, as 


dpt) = Q{t)apt ), (3.6) 

the three-dimensional metric hij in Eq. (13.21) reduces to 
hij = d?pt) (1 -l- 2CP)5ij -\- “^E^ij + Xij , (3-7) 


where 


C = C , E = E, ^ij =^ij. (3.8) 

In what follows we use an over-hat for quantities in the 
transformed frame. Erom Eq. (13.81) . the perturbations p, 
E, and Xij are invariant under the disformal transforma¬ 
tion dnj M- From Eq. (13.21) and the relation N’' = TV* 
we also obtain 

SN=isN, ip = n^ip, (3.9) 


III. DISFORMAL TRANSFORMATIONS 

In this section we discuss how background/perturbed 
quantities and the action (EH) are mapped under the 
disformal transformation EH- 


A. Transformation of background and perturbed 
quantities 


where an over-bar represents background quantities, and 


+ TX 

^ “ 02 - x^r^x ■ 


(3.10) 


Since (p = (p{t) and X = —(p^{t)/N^ in the unitary gauge, 
the quantities F and X in Eq. (13.101) contain the infor¬ 
mation of perturbations through the lapse function TV. 


B. Transformation of the action 


In the unitary gauge (EH, the line element ds^ = 
g^vdx^dx'^ in the transformed frame reads [H, HI 

dP = -N^dt^ + h^j{dx^ + N^dt){dx^ + N^dt ), (3.1) 

where 


The disformal transformation of the action Sg = 
f d‘^Xy/—g L, where L is the Lagrangia n (12.41) of GLPV 
theories, was already discussed in Refs. [d^HJ. First of 
all, the volume element y/—g transforms as 


hy = O^hy , (3.2) 


\/^= 


TV = TVa 


(3.11) 
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In the unitary gauge, the intrinsic and extrinsic curva¬ 
tures obey the transformation laws 

TZij = TZij , Kij = I ( 3 - 12 ) 


where 


a 


n 


(3.13) 


The action in the transformed frame Sg = J d^XyZ—gL 
preserves the structure of original GLPV action, such 
that 

L = A2{N ,t) + A^{N ,t)k 

+A4{n, t){k^ -s) + B4 {n, t)n 

+A5{N,t)k3 + B5{N,t) (u-^kli^ , (3.14) 

where 

~ If. 3w 6a;^ \ 


^4 = ^(^ 4 -^ 415 ], 




is = ^^s, 

• 


(3.16) 

(3.17) 

(3.18) 

(3.19) 

(3.20) 


The matter action in the transformed frame is given 
by Sm = f d*xy/—g Lm, where 




(3.21) 


where Lm depends on the JF metric g^n and the matter 
field By expressing in terms of the metric g^,^ 
in the transformed frame from Eq. dm, it contains the 
contribution of (p and its derivative. Hence the scalar field 
(j) is (kinetically) coupled to matter in the transformed 
frame. 

From Eq. (IQ) the transformation law of is given 

by 








(3.22) 


^9pi' 

On using Eqs. (11.11) and (I3.11L it follows that [s^ 

TP.^ = - T^xdf^cpd''cpd.ycpdp(P) . (3.23) 

The transformed metric with upper indices is given by 

r'' = ^ ( 9 ^'' - . (3.24) 


From Eqs. p.23l) and (I3.24E the mixed energy- 
momentum tensor obeys the transformation law: 


= TJ O a 


[ <5^ - 






djpd'^ (j)d\pd^ (j) 


(3.25) 


For the choice of the unitary gauge (IQ) . we obtain the 
relations 


rpO _ 

J- n — -in 




TZ = Tf , T* = r; . (3.26) 


We decompose the energy-momentum tensor into the 
background and perturbed parts, as Tq = —p — 5p, 
Tf = di6q, and Tf = [P SP)Sj, where Tf is a per¬ 
turbed quantity itself. The background energy density p 
and the pressure P are subject to the transformations 


^ 03^’ 


P = 




For the linear perturbations, we have 

SP ={6P - pP SN) 


where 


p = 


1 da 


N=N 


dp 

'dN 


N=N 


The quantity p is related to a and /3, as 


1 




(3.27) 

(3.28) 

(3.29) 

(3.30) 

(3.31) 

(3.32) 


In summary, the action in the transformed frame reads 

s = J (kx^f^L{N,k,s,n,U]t) 

+ j d'^x^/^Lm{gp,y{(|),dfJ,(j)),'^>m) , (3.33) 


where L and Lm are given, respectively, by Eqs. (13.141) 
and (j3.21|) . 


IV. EQUATIONS OF MOTION IN THE JF AND 
THE TRANSFORMED FRAME 


In this section we present linear perturbation equations 
of motion on the flat FLRW background for the theo¬ 
ries described by the action (1^ . We then study how 
they are transformed under the disformal transformation 

(ED. 
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A. Equations of motion in the JF 

The equations of motion in the JF were already derived 
in Refs. (4^ [s^ - fssj l for the action (12.61) . We consider the 
background line element = —N'^dt^ + a^{t)5ijdx'^dx^ 
without setting iV = 1. Defining the Hubble parameter 


From the variational principle SSg + SSm = 0, we obtain 
the background equations of motion 


L + NL^n - 377JF = p , 

(4.7) 

- .F 

L-^- 3HT = -P. 

N 

(4.8) 


H=-^, 

Na ’ 


(4.1) 


the background values of extrinsic and intrinsic curva¬ 
tures are given, respectively, by 


Kfiu = , TZgv = 0, (4.2) 

and hence K = 3iJ, S = and It = U = 0. 

The background and perturbation equations of motion 
follow from first-order and second-order Lagrangians, re¬ 
spectively, derived by expanding the action (12.61) up to 
quadratic order in perturbations. The perturbations of 
N, K, S are given by 6 N = N — N, 5K = K — 3iJ, 
and SS = 2H6K + SKj^SK’^. We write the intrin¬ 
sic curvature as TZ = SiTZ + where 5iTZ and 52 'TZ 
are first-order and second-order perturbations, respec¬ 
tively. For the perturbation W, we have the relation 
/ d'^Xy/^\{t)lA = / d'^Xy/^[\{t)TZK/2 + \{t)TZH2N)] 
up to a boundary term, where \{t) is an arbitrary func¬ 
tion with respect to t |68l |. 

We consider the perturbed metric (13.41) with the choice 
of unitary gauge 


<5^ = 0, E = 0, (4.3) 

under which the temporal and spatial coordinate trans¬ 
formation vectors are fixed, respectively. 


1. Background equations 


Expanding the gravitational action Sg = J d^x^—gL 
up to first order in the scalar perturbations, it follows 
that [13 


5S„ = / d^a 


{L + NL^n - iHF) SN 


+3a^N 3i/J'j Sa 


, (4.4) 


where 


E=L,k + 2HL,s, (4.5) 

and we dropped a boundary term irrelevant to the dy¬ 
namics. Here and in the following, the coefficients in front 
of perturbed quantities [such as those in front of SN and 
6a in Eq. (14.41) ] should be evaluated on the background. 

Variation of the matter energy-momentum tensor 
SSm = f d'^Xy/—gTf^’^6gfi,^/2 reads 


SSm = 


J d^x{-a^pSN+ 3a^NPSa) . (4.6) 


Since the matter component is not directly coupled to 
the field </> in the JE, it obeys the standard continuity 
equation 


^+3H{p + P)=0. 


(4.9) 


2. Perturbation equations 


Expanding the action (12.61) with the Lagrangian (12.41) 
up to second order in scalar perturbations and taking 
the variation with respect to SN, and the field 

perturbation Sx associated with the matter Lagrangian 
P{Y), the resulting perturbation equations of motion in 
the presence of a barotropic perfect fluid are given, re¬ 
spectively, by 


{2NL^n + N^L^nn - GHNW + UL^sH^) 


SN 

w 




W-4{NV + £)^ = Sp, (4.10) 


SN 4Ls- 

^ {a^Ny)+4iNV + £)^ + ^d\ 


(4.11) 


a^N dt 


PN a? 


SN 

-3{p + P)^=35P, 

^X+3H{Sp + SP) 

, , -dHq 


(4.12) 


(4.13) 


where 


T’ = — 


2N^ 


HL 


,NU ) 




W = L^kn + 2P[L^sN + 

^,_4L^sd^ip 
y = -5-3dg. 


4HL. 


N 


m a? 

The momentum perturbation Sq obeys 


(4.14) 

(4.15) 

(4.16) 

(4.17) 


^±{NSq)+3HNSq = -{p + P)^-^-5P. (4.18) 
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Substituting Eq. (j4.17l) into Eq. (14.121) and using 
Eq. (I4.18L it follows that 


N 


SN 


L^s 


'f 

iV2 


+ {NV + £)^+SC = 0. 


(4.19) 


where we set the integration constant 0. We define the 
gauge-invariant Bardeen potentials [nl 


^ 6 N 1 d fil; 

4/ = —^ ^—— I 
N N dt\N 


<^ = C + H 


t- 

N 


(4.20) 


and the anisotropy parameter 

d> 



(4.21) 


The effective gravitational potential associated with the 
deviation of light r ays in CMB and weak leasing obser¬ 
vations is given by [7^ 


^•efT = i(4'-4') = i(l + ??)4'. (4.22) 


One can write Eq. (14.1911 in the form 


4'-!-$ 


L,s 

NHL^S 


(C-4>) 


(4.23) 




N 


where = N^£/L^s is the tensor propagation speed 
squared discussed later in Sec. El and 


NV- 


£ ^ cl ^ NV 

““ “ L,s ~ ^ L,s 


—. (4.24) 


The parameter an characterizes the deviation from Horn- 
deski theories. Provided that one of the conditions 
L,s 7 ^ Oj ^ iV^, and Oh 7 ^ 0 is satisfied, the anisotropic 
stress does not generally vanish {rj ^ 1). 

We define the gauge-invariant matter density contrast, 
as 


<5™ = (5 - , (4.25) 


matter perturbations through Eq. (j4.26l) . The evolution 
of 41 and GeS is known by solving the coupled Eqs. (14.101) - 
(14.1311 with Eq. (14.1811 . 

On using Eqs. (14.101) and (14.111) . the second-order ac¬ 
tion of scalar perturbations can be expressed in terms 
of C, Sx, and its derivatives. Assuming that the matter 
sector does not correspond to a ghost mode, the scalar 
ghost is absent under the condition [H, 


2L,s{‘iL^sWs+3N^W^) 

7V3yp2 


> 0 , 


(4.28) 


where Ws = 2 NL,n + N'^L^nn - 6HNW + UL^sH"^- 
In GLPV theories there is a mixing between the scalar 
propagation speed Cg and the matter sound speed Cm- 
For non-relativistic matter characterized by P = -|-0 and 
5P = 4-0, we have^ = 4-0 in the small-scale limit, while 
cl is given by [H, [ssi [s^ 


c? = 


2N 

Qs 


where 


jii 4 :L^^p 


^ ^ 4Ls(JV1^ + £) 


NW 


(4.29) 

(4.30) 


This shows that the deviation from Horndeski theories 
(an 7 ^ 0) modifies the scalar sound speed. We require 
Cg > 0 to avoid Laplacian instabilities. 


B. Equations of motion in the transformed frame 


In the transformed frame described by the action 
(13.331) . we also derive the background and perturbation 
equations of motion. The cosmic time t in the trans¬ 
formed frame is related to t in the JF, as 


i = 



(4.31) 


1. Background equations 


where S = Sp/p and Vm = NHSq/p. Taking the time 
derivative of Eq. (j4.13ll and using Eq. (I4.18|l in Fourier 
space, we obtain 

P 

6 m + 2HSm + ^4/ = -3B-6HB, (4.26) 

where B = (^ + Vm, and fc is a coming wave number. We 
define the effective gravitational coupling Ges, as 

= -AnGespSm ■ (4.27) 

The gravitational potential 4' contains the information 
of gravitational coupling with the scalar field (j). The 
modihed gravitational interaction affects the evolution of 


Following the same procedure as that in 
obtain the background equations of motion 

the JF, we 

L + - 3Hf = p , 

(4.32) 

L - i - 3HP = -P , 

N 

(4.33) 

where T = L p. + 2HL ^ and 


^ d 1 dd 

~ Nd a di ■ 

(4.34) 


The relations of the quantities N, H, iF,L,L ^ with those 
in the JF are given in Appendix A (see also Ref. H). 
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Using these relations as well as the background equations 
(I4.7I) - (I4.8|) in the JF and the correspondence (13.271) . we 
can also derive Eqs. (14.32^ - 114.331) . 

The continuity equation (gH) in the JF can be ex¬ 
pressed in terms of p and P on account of Eq. (13.271) . 
Then, the corresponding equation in the transformed 
frame is given by 

§+3i7(p + P) = Q, (4.35) 

at 

where 




H{\ — afi) + afi— 

7VJ 


P. 


(4.36) 


The rhs of Eq. (14.351) describes the coupling between mat¬ 
ter and the scalar degree of freedom. 


2. Perturbation equations 


In the transformed frame, the perturbation equations 
following from the variations with respect to 5N, i9^'0, 
and are given, respectively, by 

(2NL j^ + - 6HNW + 12L ^ 

^3^ - NW-4{NT> + £)^ =Sp, (4.37) 


dt 

N N dt 

cP di 

^ 6N 

-3(p + P)V =3JP, 
N 


44. (d^Ny) + + f)4^ + 

dt\ > n2]\T 


a?N 




(4.38) 


(4.39) 


where I),S,W,y are defined in the same way as 
Eqs. (I4.14I) - (I4.17I) with the over-hat for background 
quantities. We can also derive Eqs. (14.3711 - (14.391) from 
Eqs. (I4.10I) - (I4.12I) in the JE by using the relations (13.8L 
(13.91) . (lAll) and (IA2F 

The equation of motion for the momentum perturba¬ 
tion Sq can be derived by employing Eqs. (I3.9|l . (I3.27L 
( 12211 ), (ESI, (OSl) - and dSH) , as 


4 {N6\) + SHNSq = -{p + P)^-SP. 
dt TV 


(4.40) 


This is of the same form as Eq. (14.181) in the JF. 
Similarly, the density perturbation Sp obeys 

P di 

+mp [H-^][5P + pPpsn'^ 


d 

di 


[Sp — i'pSn'J 


ZH- 


(^p — vpdN'^ 


LO 

N 


-k 




d^N. 




the form of which is modified relative to Eq. (14.131) in the 
JF. This modification comes from an explicit coupling 
between matter and the scalar degree of freedom in the 
transformed frame. 

From Eqs. (I4.39|) and (14.401) it follows that 



+ ~ 


N dt ’^J 


i. 

N 


+ {M + £)^+SC = 0 . 

N 


d^,s dp) 

N di 


(4.42) 


We introduce the gauge-invariant Bardeen potentials in 
the transformed frame, as 


N dt \nJ 


4 = C -k ni . (4.43) 

N 


From Eq. (14.421) we obtain the relation similar to 
Eq. (j4.23|l with additional over-hats to each quantity. As 
we see in Sec. El it is possible to find a metric frame in 
which some of the terms generating the anisotropic stress 
vanish. 

We also introduce the gauge-invariant matter density 
contrast, as 


5^ = 6-ZVm, (4.44) 

where 6 = dp/p and Vm = NHSq/p. From Eqs. (j4.40l) 
and (14.411) we can derive the second-order equation for 
5m analogous to Eq. (14.261) in the JE. If we transform to 
the EF, the effective gravitational coupling with matter 
becomes particularly transparent. We shall address this 
issue in Sec. El 


V. EINSTEIN FRAME 

We dehne the EF in which the second-order action of 
tensor perturbations "fij is of the same form as that in 
GR gj. In the following we discuss the transformation 
of the action in the JF frame to that in the EF. 


A. Transformation to the EF 

Expanding the action (12. 6 |) with the Lagrangian (12.41) 
in terms of tensor perturbations 7 ^, the resulting 
quadratic action of 7 ^ in the JE is given by 

= J d'^xa^qt5^^5^^ (iijiki - , (5.1) 

where 



4^ 




a 


(5.2) 



























































In Eq. (EH) the quantities qt and Cj should be evaluated 
on the background, such that the kinetic term X appear¬ 
ing in L^s = —Ai — and £ = + B^/{2N) cor¬ 

responds to the time derivative X{t) = —(jP‘{t)/{2N‘^). 
We require the conditions qt > 0 and > 0 to avoid 
ghosts and Laplacian instabilities. In GR we have L = 
Ml,R/2 = -{M^J2){K^ -S) + (MpV2)7^, in which case 
qt = XI^^/{8N) and (where Mpi is the reduced 

Planck mass). 

Under the disformal transformation (EH the tensor 
perturbation is invariant, see Eq. (13.81) . Hence the 
second-order tensor action in the transformed frame 
reads 


= 


d'^xa^qt5^^5^'- 




iijiki - ■^djijd'jki , (5.3) 


where qt = L ^/{AN) and = N^£/L g. Comparing 
Eq. ()5.3I) with Eq. (15.ip and using the relation (j3.6l) . it 
follows that 


. _ 1 

qt - -^qt , 


= 


(5.4) 


The EF corresponds to a frame in which both qt and 
are of the same forms as those in GR, i.e., qt = M'^J[8N) 

and = 7V^. The tensor action EH in the JF can be 
transformed to that in the EF for the choice 


= 


8gtCt 

Ml 


r = 


SgtCt 

Ml 


A 

N2 


- 1 


1 

X 


(5.5) 


The quantities qt and Cj in the action (15.11) depend on 
the time t alone. Then, the factor T in Eq. (15.5p has 
the dependence r(()), X) = 'y{(j))/X in the unitary gauge, 
where = {8qtCt/Ml){c^/N^ - 1 ). 

For the choices (15.51) the terms a and /3 in Eqs. (13.31) 
and ()3.10l) are given by 


1 r.Ct 

a = — = . 

P N 


(5.6) 


Since both a and P are functions of t, we have ^ = 0 = u 
from Eq. (j3.3ip . Then, the coupling (14.361) reduces to 


Q = -^f/3-3P) 


N 


NCt 


where 


^ _ ct 


(5.7) 


(5.8) 


B. Background equations in the EF 


The choice (15.51) corresponds to the conditions 

Ml 

L g = -A4 - SHA^ = , 


£ = B4 


1 dB^ Ml 

2 dt 2 


(5.9) 

(5.10) 


Using the relation (EH, the background Eqs. (I4.32p and 
(14.331) in the EF can be written in the following forms: 

8MlH^ = pj,E + P, (5.11) 

dH 

—2Mpi-^ = Pde + A)e + P + P, (5.12) 


where 


Pde = —A2 — 6 H H5 
-N (4,^ + 
Pde = ^2 + 


5,N I ’ 


(5.13) 


_ ( ^ _ UH^A^ - 6P2^ ) . (5.14) 
\ dt dt dt ^ ' 


From Eqs. (13.271) and (15.61) . the matter equation of state 
w = P/p is invariant under the transformation to the 
EF, i.e.. P/p = P/p. The energy density pde and the 
pressure Pde obey the equation of motion 


dpDE 

di 


P8H 


(^Pde + Aje) = —<3, (5.15) 


where Q is given by Eq. EH- Comparing Eq. (14.351) 
with Eq. (I5.15P , the scalar field and matter interact with 
each other in the EF. 

The background equations of motion in the JF do not 
contain the terms P 4 and Pg. The theories with same 
values of H. 2 , 3 , 4 ,g but with different P 4 _g cannot be distin¬ 
guished from each other at the background level l5i,[73. 
In other words, two theories with different values of 
lead to the same background dynamics for given H. 2 , 3 , 4 ,g. 

This implies that the coupling —Ct/{NCt)p appearing 
in Eq. EH does not essentially modify the background 
physics even for the theories in which Ct varies in time. 
In Sec. El we shall confirm this property for a concrete 
dark energy model. 


C. Perturbations in the EF 


If Ct = N, then T = 0 and Q = —{u)/N){p — 3P). This 
case corresponds to the well-known conformal transfor¬ 
mation arising e.g., in BD theory [^ . For radiation 
(p = 3P) the coupling Q vanishes, but for non-relativistic 
matter (P = 0), we have that Q = —{uj/N)p. If Ct varies 
in time, the last term on the rhs of Eq. EH does not 
vanish even for radiation. 


Substituting the relations L g = £ = Ml/2 into 
Eqs. (14.371) - (14.421) . we obtain the perturbation equations 
of motion in the EF. From Eq. (I4.42p the gauge-invariant 
Bardeen potentials obey the relation 


7 A - 

4'-k 4* =-an — , 
N 


(5.16) 
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where dn = 2 N'D/M^^ is the parameter characterizing 

the departure from Horndeski theories. Since L ^ = 0 

and Ct = the first and second terms present on the 
rhs of Eq. (14.231) in the JF vanish in the EF. 

The full GLPV action cannot be mapped to the full 
Horndeski action under the disformal transformation 
[ 4 ^ , so the parameter dn in Eq. (I5.16|) does not gen¬ 

erally vanish. It is, however, possible to transform part of 
the GLPV action to the action belonging to Horndeski 
theories, in which case dn = 0 and hence there is no 
anisotropic stress in the EF. 

Let us consider perturbations of non-relativistic mat¬ 
ter characterized by P = 0 and SP = 0. In the EF, 
Eqs. (14.401) and (14.411) reduce, respectively, to 


I dV^ I 1 d/3 1 dH\ ^ _ SN 
H di dt m di ) N ’ 

d6 dC k'^ 1 

^ + 3-^ + = 0 , 

dt dt a? jy 0 ? H 


(5.17) 

(5.18) 


where we used the background Eq. (14.351) . Taking the 
t derivative of Eq. (I5.18|) and employing Eq. (15.171) . the 
matter density contrast (14.441) obeys 


d'^dm 

dP 


= -3 


2H- 


1 dfl 
di 


Ct di 


dSjri 

di 




d'^B 

dP 


1 dit 1 dCt \ dB 


hi- 


where B = (^ + Vm, and 




/ 1 dfl 
di 


Ct di ) 


at. 

N 


- 1 


(5.20) 

N 


The effective potential characterizes the strength of 
gravitational coupling with matter. In the EF, it is clear 
that 'kg is expressed in terms of the sum of the gravita¬ 
tional potential d)* and contributions from the variations 
of H and Ct as well as the difference of Ct from 1. For 
the theories with Ct = 1, which is the case for BD the¬ 
ories, the variation of the conformal factor H gives rise 
to the modification to 'k. The deviation of Ct from 1 
and the variation of Ct occur for the theories studied in 
Refs. [ 5 ^ [ 15 , in which case the gravitational interaction 
is modified as well. 

On using the correspondence SN/N = SN/N and 
t/}/N = Ctt/}/{NCt), the gravitational potential 4' can be 
expressed by using 4* in the JF. Then, we obtain the 
simple relation 

^g = J- (5.21) 


This shows that the effective potential 4/g is directly 
related to 4' appearing in the matter perturbation 
Eq. (14.261) in the JF. Once we find the evolution of 'kg 
in the EF, the potential 4' and the effective gravitational 
coupling GeS in the JF are known accordingly. 


D. Model belonging to Horndeski theories in the 

EF 


One example of realizing du = 0 in the EF is the model 
described by the JF Lagrangian 

L = A 2 (TV, t) + H 3 (TV, t)K + A 4 {t) {K^ - S) + B^ {t)n. 

(5.22) 

In this case the tensor propagation speed squared di¬ 
vided by TV^ reads 


(5.23) 

/I 4 

In Horndeski theories we have —^44 = P 4 from the first 
condition of Eq. ( 1 ^ . but in GLPV theories there is no 
such restriction and hence C^ generally differs from 1 . 
Since P = 0 and — 1 for the model (15.221) . it 

follows that 


4 / -p $ 


H 4 

TVPH4 


(C - $) - (c^ - 1 ) 



(5.24) 


from Eq. (14.231) . For A 4 depending on t and for C^ dif¬ 
ferent from 1, the anisotropic stress is present in the JF. 
Under the disformal transformation with the factors 




2v^—H4P4 


2-^/— ^4^4 



the Lagrangian (15.221) is transformed to 


(5.25) 

(5.26) 


L = A2(N,t) + A3(N,t)3C + ^ +TZ') , 

(^.•27) 

where 242 = {A 2 — ‘iujA 3 /N+%uj'^A 4 /N'^)/{CPa) and Ha = 
[A^ — 'iijjAi/N)/^^ with a = (•\/2/TWpi)(— 54 / 414 )^/^ and 
Lo = {A 4 IA 4 -f B^lBPjjA. The last term of Eq. (j5.27l) 
corresponds to the Einstein-Hilbert term TV/pjR/2, where 

R is the four-dimensional Ricci scalar. Since dn = 0 for 
the Lagrangian (j5.27|) . it follows that 


^-bd = 0. (5.28) 

Thus, for the Lagrangian (15.221) . the disformal transfor¬ 
mation allows one to de-mix the gravitational potentials 
in the EF, such that the anisotropy parameter fj = —4>/'k 
is equivalent to 1. 

While the gravitational potentials are de-mixed in the 
EF, the matter perturbation Sm is subject to gravita¬ 
tional mixing described by the effective potential (15.201) 
mediated by the scalar field (/)■ Multiplying the term 
3NH for Eq. (14.381) and taking the sum with Eq. (I4.37L 
one can relate dm with metric perturbations. Let us 
employ the sub-horizon approximation under which 
the dominant contributions to the Ihs of Eqs. (j4.37l) 
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and (I4.38|) are those involving and d'^((a^. In 

the EF the terms £ and >V are given, respectively, by 
£ = Mpj/2 and VV = 2HM^^/N + Provided that 

2i?Mpj/7V ^ 1^3 jvl’ obtain the Poisson equation 



1 


p5j^ 


(5.29) 


On using Eq. (15.281) and introducing the gravitational 
constant as G = ( 87 rMp[)“^, it follows that 


—4' ~ -AirGpSm ■ (5.30) 

This shows that the gravitational coupling associated 
with is simply given by G under the condition 
2HM^,/N » 14^^10 Hence, the modified gravitational 
interaction from GR arises from the terms on the rhs of 
Eq. (15.201) other than T. 

In Sec. EH we consider a concrete model belonging to 
the Lagrangian (15.221) and study the correspondence of 
physical quantities between the JF and the EF in detail. 


VI. CONCRETE MODEL 

We study a dark energy model described by the action 
where the Lagrangian L is given by Eq. (15.221) . We 
consider the following functions 

A2 = -^e{^)X-Vi(^), Ha =-Mp2 , 

Ai =, 54 = iMp2E2((/)), (6.1) 

where e((/)), V{cj)), E’i(^), and F2{4>) are functions that 
depend on cj), i.e., on t in the unitary gauge. The de¬ 
pendence of A 2 and H 3 on N arises in the kinetic term 
X = — To accommodate BD theory as well 

as theories recently proposed in Refs. [^, [^, we choose 
the functions 

Fi((^) = 6-291, F2(((>) = 46-29=-^/^^' , 

e((/)) = (l- 6 g?)Fi(<^), ( 6 . 2 ) 

where qi,q 2 ,cti are positive constants. We assume that 
(/i, 92 1 for the compatibility with observations [^, 

It^ . In Horndeski theories, the first condition of Eq. (1231) 
demands that F2{4>) is equivalent to Ei(</>). The original 
BD theory without the field potential corresponds to the 
case V{(j)) = 0 and F 2 (<(') = Fi {(p) with the BD parameter 
WBD = (1 - 6qi)/Aqf) [ 13 . 


^ There are some models Hire Irinetic braidings [T3l in which the 
dependence of A 3 on N modifies the gravitational coupling; see 
Refs. [Ii,[23. 


Let us first consider theoretical consistent conditions 
in the JF. In the following we set the background value 
of the lapse N to be 1. From Eqs. (15.21) and (I4.28|) the 
conditions for avoiding tensor and scalar ghosts are given, 
respectively, by 


qt = > 0 , 

(6.3) 

M2^4Fi ^ ^ 

9s = -:— > 0, 

2{Mp,H-qi<PY 

(6.4) 


which are satisfied for Fi > 0. The tensor and scalar 
propagation speed squares are given, respectively, by 

C? = = 0^2(91^ (g 5) 

2 2 , (l-c 2 )Dm 2 {qi - q2)ci{\/6 - 6qiXi) 

+ 2x2 +-3^^-’ 

( 6 . 6 ) 


where 


Xi = 


y/6HMpi 


U 777 , — 




(6.7) 


Under the no-ghost condition Fi > 0, the condition 
(HID is satisfied for ^2 > 0. For the theories with gi = 92 , 
c2 is constant {cf = c^J. Since Cg = c2-|-(l —c2)Om/(2a;2) 
in this case, Cg is positive for 0 < < 1 , while Cg can 

be negative for > 1. In Ref. [ 5 ^ the authors studied 
the cosmology for the specific case with 91 = 92 = 0. If 
9 i A then varies in time. The variation of gives 
rise to a contribution to Cg, i.e., the last term on the rhs 
of Eq. (16.61) . The cosmo logy with 91 = 0 and 92 ^ 0 was 
recently studied in Ref. j56l[ as a model of realizing weak 
gravity on scales relevant to large-scale structures. 

For dark energy models in which the ratio £lrnlx\ de¬ 
creases with time, Cg grows to be very much larger than 1 
as we go back to the past. This behavior can be avoided 
for the scaling model characterized by the potential [t^ 

V{p,) = -b , (6.8) 

where Vi, V 2 , Ai, A 2 are positive constants. Provided the 
first potential on the rhs of Eq. () 6 . 8 I) dominates over the 
second one, the scaling solution with the constant ratio 
is realized during radiation and matter eras 
The solution exits from the scaling regime to the epoch 
of cosmic acceleration due to the existence of the second 
potential. We shall consider the situation in which the 
slopes Ai and A 2 are in the range 

Ai>10, A2<1, (6.9) 

for consistency with the big-bang nucleosynthesis [t^ 
and the late-time cosmic acceleration [l[ . There are seven 
free parameters ( 91 , 92 , Cti, Ui, V 2 , Ai, A 2 ) in our model. 
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A. Transformation to the Einstein frame 

For the theories given by the functions (EH), the two 
factors ()5.25l) and (I5.26P transforming to the EF are 
given, respectively, by 

= Ei(0)ct(<^), (6.10) 

mX) = y^F^F^) ( 1 ^ ■ 0 Y 


The disformal transformation to the EF corresponds 
to the change of tensor propagation speed squared = 
—B 4 /A 4 to Cj = 7V^, so the quantity q 2 arises in the EF. 
Nevertheless we are dealing with the same physics in the 
two frames, so any physical condition (such as the stabil¬ 
ity of fixed points) should not be subject to change. In 
what follows we shall study the correspondence of back¬ 
ground quantities between the EF and the JF. 


We also have 


a = ^ = \J■ (6.12) 


Then the action in the EF is given by Eq. (13.331) with the 
Lagrangian (I5.27L where 


^2 — 
^3 = 


2iV2 

N 


1 - 2 ( 91 - 92 )" -!"(</>), (6.13) 


(91 - 92 ): 


and the EE frame potential 


E(0) = 


!/(</.) 




,(gi-l-392)</i/Mpi 


For the JF potential (16.81) it follows that 

(>(,),) = -f , 

where Vi = V 2 = and 


(6.14) 


(6.15) 


(6.16) 


1. Einstein frame 

The background equations of motion in the EF are 
given by Eqs. (I5.11I) - (I5.12|) . where 


1 f d(j) 


= 2 [Ti 


1-^(9i-92)" 


!>(</.) 


-1-3(gi - q2)MpiH 


def) 

di 


( 6 . 21 ) 






i-jta-®)" -r(« 

d'^cj) 


( 6 . 22 ) 


The matter fluid and the scalar field cj) obey Eqs. (I4.35() 
and (j5.15p . respectively, with 


Q = 


9i + 92 7 . 9i - 92 . d(j) 


2M, 


pi 


(;i-3P) 


di 


Mp\ di 


(6.23) 


Hi — Xi — qi — 3(72 , H 2 — X 2 — qi — 3^2 • (6-17) 

Eor the theories with qi = 52 (be., = constant) the 

term A 3 vanishes, in which case the Lagrangian (15.271) 
corresponds to that of the canonical scalar field 4 > coupled 
to matter. 


B. Background dynamics 

Erom Eqs. (I4.7I) - (I4.9I) the background equations of mo¬ 
tion in the JF are given by 

3M^^H^Fi = i^^ + V- 3M^,H^Fi^^ + p, (6.18) 

-2M^iFiH = e<P^ + M^iiFi - HFi) + p + P, (6.19) 
p-f 377(p-f P) = 0 . (6.20) 

These equations depend on the function Fi[<f>) but not 
on F 2 { 4 >), so the quantities 52 and are irrelevant to 
the background dynamics. This means that the theories 
with same qi and different (72 (i.e., same A 4 and differ¬ 
ent B 4 ) cannot be distinguished from each other for the 
background cosmology in the JF [s^. 


The first term on the rhs of Eq. ()6.23p arises for the 
standard coupled dark energy scenario characterized by 
91 = 92 and Cjj = 1 [H^. The second term on the rhs of 
Eq. (I6.23|) does not vanish for the theories with qi ^ 52 
(i.e. time-varying c^). 

To discuss the background dynamics in the EF, it is 
convenient to introduce the dimensionless quantities 


1 


Xi = 


d(j) 


Ude = 


\/&Mp\H dt 
Pde 


X2 = 




y/iMp\H 


3M2jiJ2 


Cl = ^ 


3Mp2 772 ’ 


P 

w = — . 
P 


H = 


MpiV^ 

V 


(6.24) 

H 


where a prime represents a derivative with respect to 
Af = J Hdt. Recall that the matter equation of state 
w is invariant under the disformal transformation to the 
EF. In the following we assume that w is constant. From 
Eq. (15.111) there is the relation IIde + Clm = 1- 

The variables xi and X 2 obey the differential equations 
















































12 


/fi 

x'l = —[2qi(3w - 1) - 2y6{l + gi(gi - g2)(3w - l)}ii + gi{3(gi - 52 )^ - 2}(3ii; - l)x\ 
+2{3g2 + gi(2 + 3 w)}x2] + xi^h , 

X2 = - ^IJ'XiX2 + X2eH , 

3 3 3 

iH = ^{1 + W - xl) - -qiiqi - g2)(3w - 1) + ■^VQiqi - g 2 ){l - w + gi(gi - g2)(3w - l)}ii 
-^{3(gi - 92 )^ - 2}{1 - w + gi(gi - g2)(3w - l)}xj - - (gi - g2){gi(2 + 3w) - 3g2 


(6.25) 

(6.26) 


fi}]xl , (6.27) 


where in is related to the effective equation of state WeS of the system, as WeS = — 1 + llujZ. The field density 
parameter and the equation of state are given, respectively, by 


floE = 1 - flm = ^ 2V6(gi - 52) + {2 - 3 (gi - 92)^}®! 


2 L 


^2 5 


WbE = 


Pde _ 3[2 - 3 (gi - 92)^]:^? - 6*2 - 2^6(91 - 92 )(i'i - ch^i) 


Pde 


3[2 - 3(91 - 92)^]if + 6 x 2 + 6v^(9i - 92)^1 


(6.28) 

(6.29) 


If n is constant, which is the case for the exponential 
potential V{cj)) = there are five fixed points 

characterized by constant xi and X 2 - They are summa¬ 


rized in Table m 

For the fixed point (a), the field density parameter 
reads 


9 i(3w - 1 )[9 i{4 - 3(91 - 92)^} - 692 -f 3{3gi(9i - 92)^ -b 292}^] 
3[1 - w -I- 91(91 - 92 )( 3 u; - 1)]2 


(6.30) 


so that both Ude and xi vanish for radiation (w = 1/3). 
If 92 = 9 i, it follows that xi = y/&qi{3w — 1)/[3(1 — w)], 
WeS = [29i (1 — 3w)^ -I- 3?u(I — w)]/[3(I — w)], and I/de = 
29 ^( 3 w — l)^/[3(r(;— 1)^] for the point (a). When w = 0, 
this corresponds to the (/-matter-dominated era (^MDE) 
characterized by zieff = ^de = 2qf/3. Provided that 
w ^ 1/3, the point (a) is a kind of scaling solution with 
a constant ratio flm/^DE- 

Since the effective equations of state WeS for the points 


(bl) and (b2) are I, they are neither relevant to radia¬ 
tion/matter eras nor the late-time cosmic acceleration. 

The point (c) is the scalar-field dominated point 
(Ude = I) relevant to dark energy. When 92 = 91 we 
have x\ = X 2 = ^1 — Ai^/ 6 , and Weff = —I-l-p.^/3, 

so the cosmic acceleration occurs for y? <2. For 92 ^ 91 , 
z&efi is close to — 1 provided that 91 , 92 , /r are smaller than 
the order of 1 . 

The point (d) corresponds to the scaling solution with 
the field density parameter 


ql(9w'^ — 30w — 23) -|- 291 ( 8/1 -I- 392 {9 -|- w{‘i + 3z(;)}] -I- 3{w -I- 1)[4 — 492 /t -I- Sq^iw — 3)] 

[ 2/1 - 91 - 1 - 392 - 3w{qi + 92 )]^ 


(6.31) 


When 92 = 9 i, we have that Ode = [3(1 +w)— qi{3w — 
l)(/i -I- 9 i - 3qiw)]/{^ + 9 i - 3qiw)‘^ and zieff = [/izc -b 
qi{3w — !)]/[/! — qi(3w — 1)]. In this case the radiation 
scaling solution corresponds to Ode = and zieff = 
1/3, whereas the matter scaling solution is characterized 
by Ode = {3+qin+ql)/{^j,+qi)'^ and zieff = -qi/{fi+qi). 
If the fixed point (d) is stable, the solutions approach it 
during the radiation and matter eras. 

The stability of the above fixed points can be analyzed 


by considering small perturbations 5xi and 8 x 2 about 
each of them. The linearized version of Eqs. (I6.25|) and 
(16.261) can be written in the form 



where A4 is a 2 x 2 matrix. If the two eigenvalues ki _2 of 
M. are negative (or imaginary with negative real parts), 
then the corresponding fixed point is stable. 
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Xl 

xl 

foeff 

SJde 

(a) 

•\/6(4i(3ui—1) 

0 

3444(1 —444)+9l (3444—1) [91(610 + 1) —392] 

Ea.dODl) 

3 [ 1 - to-1-(Jl ((J1 - (J2 ) ( 3 444 -1) 1 

3(l-444) + 39l (9l-92)(3l44-l) 

(bl) 

\/6 

VE+3{qi-q2) 

0 

1 

1 

(b 2 ) 

\/6-3((?i-(72) 

0 

1 

1 

(c) 

\/6(3(ji-3(j2-44) 

2[6-(39i-392-44)^] 

6-(391-392-2/4)(39i-392-/4) 

1 

3[(3(ji-3(72-/4)(qi-(?2)-2l 

3f2-(9i-92)(39i-392-7t)F 

6-3(91-92) (391-392-44) 

(d) 

^6(1+444) 

2 [ 3 — 3 w '^— 2 qi { 3 q 2 -\- fJ ^){^ w — l )-\- 2 qf { 9 w '^-\- 3 w — 2 )] 

91+392(414—1)+(39i+244)io 

Eq. (pmi 

2/4-(Ji+ 392-3to((7i-1-92) 

[2/4-91+392-310(91+92)1^ 

2/4-91+392-3(91+92)444 


Table I. The fixed points in the EF and corresponding values of Weff and IIde for the system characterized by the autonomous 
Eqs. (I6.25I) - (I6.26I) with Eq. (16.2711 . The scaling radiation and matter points (d) are stable for /r > 10 and qi,q 2 <C 1, whereas 
the accelerated point (c) is stable for /r < 1 and qi,q 2 1. The potential (16.161) allows for the transition from the matter point 
(d) with /r ~ > 10 to the point (c) with ^ ~ /r 2 < 1 . 


In the presence of radiation (w = 1/3), the eigenvalues 
of the point (d) are given by 

.(d) M - 3gi + 3(72 ± \/64 - 15(/r - 3qi + iq 2 Y 
2(7t-gi+g2) 

(6.33) 

For /i > 10 and qi,q 2 “C 1, imaginary with neg¬ 

ative real parts, so the point (d) is a stable spiral. For 
non-relativistic matter {w = 0), the eigenvalues of the 
point (d) read 


.(d) _ +3g2)± \/D{d) 

2{2n-qi+3q2) 


(6.34) 


where Zl(d) = 9 (/i - qi + ^ 2 )^ - 24[3 - 'iql + 2qi{ii -|- 
392 )][(/^ +3 ( 72 )^ - 6(71(7/-I- 3 ^ 2 ) -3 + Qq\]. For 7 / > 10 and 
'? 1 j 92 “C 1, the eigenvalues (16.341) are again imaginary 
with negative real parts. Thus, the first potential on the 
rhs of Eq. (16.161) leads to the scaling radiation and matter 
eras driven by the fixed point (d) with 74 = 7 / 1 . 

The point (a) can be potentially relevant to radiation 
and matter eras, but one of the eigenvalues is positive, 
i.e., fSf’ = 2 for w = 1/3 and = (3 -|- 2qi^ + 651^2 — 
4 ( 7 ^)/[ 2 (l-|-( 7 ig 2 — 9 i)] for w = 0. Hence the solutions are 
attracted by the scaling solution (d) rather than the point 
(a). For 74 smaller than the order of 1, the stable scaling 
matter solution (d) with Ode < 1 does not exist [t^, in 
which case the matter era is replaced by the ^MDE (a) 
[53 ■ In our model, we do not consider this latter case 
to avoid very large values of Cg in the early cosmological 
epoch. 

After the dominance of the second potential on the rhs 
of Eq. (I6.16|) . the solutions exit from the scaling matter 
era (d) to the epoch of cosmic acceleration driven by the 
point (c). In the presence of non-relativistic matter, the 
eigenvalues of the point (c) are given by 


(c) 

(c) 


1 


6 - (m - 3 gi-b 3^2)^ 

2+(gi-92)(M-3gi+3g2) ’ ^ ^ ^ 

6 — 2(74 — 3qi + 3(72)(7i — 2qi + 3 ^ 2 ) /„ g,., 

2 + (gi-g2)(7^-3(7i+3(72) 


For 74 < 1 and qi,q 2 ^ 1, it is clear that both k) ' and 

/s fc') 

4/2 ' are negative. Hence, the solutions finally approach 
the accelerated attractor (c) with w)eff close to —1 for 
7 i 2 “C 1 (see Table U for the value of Weff)- 

In summary, for the potential (I6.16|) with 741 > 10 and 
7 f 2 ^ 1, the background cosmological sequence in the EF 
is as follows: (i) scaling radiation point (d) with 44 ; = 1/3 
and 74 = 741, — >■ (ii) scaling matter point (d) with w = 0 
and 74 = 741, —>■ (hi) accelerated point (c) with 74 = 742. 


2. Jordan frame 


The background dynamics in the JF can be known by 
using relations for physical quantities between the two 
frames. We define the dimensionless quantities 


Xi = 


VQMpiH ’ 


_ MpiV^ 

A= ^ 


X 2 = 


n 


m — 


^JWiMpiH ’ 
P 


(6.37) 


where the field density parameter is given by Hde = 
1 — On using Eqs. (13.271) . (I6.10|) . (16.121) . and (lAll) . 
we obtain the following correspondence 


Xl = {I + U:h)xi , X2 = (1 -b Wi/)X2 , 

Om = (1 + W/4)^Om , (6.38) 


where 


h y/Q '/Q{qi+q2)xi 

" 440 ==-iwitaWiF ■ 

(6.39) 

The slow-roll parameter tfj = which is associ¬ 

ated with the effective equation of state 4(;eff in the JF as 
WeS = —1 -b 2ei//3, satisfies the relation 


Cff = (l-bw//) 



3q2)xi + 


1 duJH 
Hil + con) di 
(6.40) 
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The slope A defined in Eq. (16.371) is related to the slope 
/i in the EF, as 

= X- qi-3q2. (6.41) 

Using the above correspondence, one can readily trans¬ 
late the background cosmological dynamics in the EF to 
that in the JF. In the JF, the fixed point (d) in Table U 
corresponds to 

(d) _ ^6(1 + w) 

~ 2A ’ 

^(d) ^ v^3(l - w^) + 2gi(l - 3w)[A - 3gi(l -b w)] 

V 2 X ’ ■ 

with 


(c) with A = A 2 . During this sequence, the effective 
equation of state and the field density parameter evolve 
as (i) Wes = 1/3-8qi/(3Ai), Dde = 4(l-4qf-h2qiAi)/A? 
(radiation era), —?► (ii) WeS = —2(j'i/Ai, Dde = [3(1 — 
Aql) + lqi\i\/\\ (matter era), —>■ (in) WeS = —(3 — A^ — 
20qJ -I- 9 giA2)/[3(l - Aqf + q'iA 2 )], Dde = 1 (accelerated 
era). 

C. Perturbations and matter-scalar couplings 

We consider the evolution of cosmological perturba¬ 
tions and the resulting matter-scalar coupling in the 
presence of non-relativistic matter satisfying P = 0 and 
6 P = 0. In what follows we shall focus on the case where 


Wgg = w - 2(1-I-w)^ , (6.43) 

A 

^(d) 3(l-bw)(l-4gJ)-h(7iA(7-f 3w) 

“DE ~ \2 ■ 


qi=q 2 = q, (6.49) 

under which is constant. Then, at the background 
level, the coupling (16.231) reduces to 


The fixed point (c) translates to 

(c) _ A - 4gi 

^ a/6(1 - Aq1 + qiX) ’ 

(c) _ \/l - (A - 4gi)^/6 
^2 1 - 4g2 + q,X ’ 

with 

(^c) 3 — — 20^^ -l- 9^iA 

^ 3{l-4qf + qiX) 


(6.45) 


(6.46) 

(6.47) 


The quantity q 2 disappears after the transformation from 
the EF to the JF, which reflects the fact that the factor 
F 2 is absent in the background equations (16.181) - (16.201) . 

The stability of fixed points should be independent of 
the values of 92 - Substituting Eq. (16.411) into Eqs. (j6.33l) - 
(16.361) . it follows that the numerators of the eigenvalues 
do not contain the term 52 ■ In the denominators there are 
still g 2 -dependent terms, but they are simply associated 
with the transformation of the number of e-foldings, i.e., 


dAf 

^ = 1 -—(?. + «)«- 


(6.48) 


The evolution of homogenous perturbations bxj oc 
{j = 1, 2) in the EF translates to the JF evolution pro¬ 
portional to where Kj = Kj[l — V&iqi + q 2 )xi/ 2 ]. 

On using the index kj, the <72 dependence in the denom¬ 
inators of kj vanishes identically. Provided the rhs of 
Eq. (16.481) is positive, which is the case for 51,(72 "C 1 
and |a;i| < 1, the stability conditions of fixed points are 
identical to each other in the two frames. 

The above discussion shows that, in the JF, the scaling 
radiation fixed point (d) with in = 1/3 and A = Ai is 
followed by the scaling matter point (d) with w = 0 and 
A = Ai, and then the solutions finally approach the point 


^ qp d(l) 

Mpi di 

The quantity ujh in Eq. (16.391) reads 

Fi a $ q(j) 

^ 2HFi ^ ~lip ^ ~ HMpi ■ 


(6.50) 


(6.51) 


In the EF, the gauge-invariant matter density con¬ 
trast is defined by Eq. (j4.44l) . Since 5 = 6 and Vm = 
{l+ujH)ym, 5m is not equivalent to 5m- For the perturba¬ 
tion deep inside the Hubble radius the velocity potential 
Vm is much smaller than 5, so the difference between 5m 
and 5m is small. We introduce the effective gravitational 
coupling Geff in the EF, as 

= -47TGeffp5m , (6.52) 


where d'g is given by Eq. (15.201) . For the choice (I6.49L 
'kg reduces to 


qd)' ^ 




(6.53) 


where 


X = 


H'^ 

N 


(6.54) 


Using the relation (15.211) as well as the approximation 
5m — 5m for the perturbations deep inside the Hubble 
radius, we can rewrite Eq. (I6.52p of the form (/c^/a^)'!' ~ 
—4Tr{Gee/Fi)p6m- Hence, the effective gravitational cou¬ 
pling Geff in the JF is related to Geff, as 


Geff 


Geff 

Fi 


(6.55) 
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The above discussion shows that, once '^g is known by 
solving the perturbation equations in the EF, the gravi¬ 
tational potential tj/ and the resulting matter-scalar cou¬ 
pling Geff in the JF are determined accordingly. 

Since oh = 0 in the EF, we have the relation (I5.28L 
i.e., 

$ = -4'. (6.56) 

On using the relations SN/N = SN/N, C = C; and x = 
H-ip/N = Cjx/(1 + loh) with TV = 1, the gravitational 
potentials 'k and d> in the JF are related to and $ in 
the EF, as 

Vk = d- + (c? - 1) ^ - c?-^^x, (6.57) 

^ ’ dt\H) + 

^ = -(6-58) 

\1 + UJH J 

The anisotropy parameter 77 in the JF generally differs 
from 1 due to the presence of the perturbation x- 


1. Einstein frame 


where 

Cl = [2(3ii - V 6 n)ixl - 1) 

+Clm{5xi — 2V6{tJ, + q)}]/(2xi) , (6.65) 

C 2 = 3[12ii — 4'/6fj,xf — 2xf(9 + 2q^ + 3qn — 60^) 
+ixi{2q^ + — 1 — 2g(g J- 77 ))} 

- V^Xi{4g -4/7-1- (3g -I- 5/t)Om} 

--/6(g-I-/7)Om(Om - l)]/(2ii). (6.66) 

The general solution to Eq. (|6.64p can be written in 
the form Vm = Vm'^ + Vm \ where Vm'^ is the special 
solution and Vm'^ is the homogenous solution derived by 
setting the rhs of Eq. (16.641) to be zero. For the sub¬ 
horizon perturbations satisfying K. ^ 1 , the homoge¬ 
nous solution induces the rapid oscillation of the veloc¬ 
ity potential with a frequency associated with the term 
(1 - c?)flm/C^/(2c?xf). 

For the theories with g = 0, as long as the homogenous 
solution is initially suppressed relative to the special so¬ 
lution, it was found in Ref. that the perturbations 
X, C) Vm in the JF stay nearly constant during the 
scaling matter epoch. As we confirm later numerically, 
this is also the case for g yf 0. At the scaling fixed point 
(d) with w = 0 and fi = fii, the ratio Ctm/xl is constant. 
We consider the sub-horizon perturbations satisfying the 
condition (1 — /{^c^xD » IC 2 I. Then, the spe¬ 

cial solution to Eq. (I6.64|) is given by 


Let us study the evolution of perturbations in the EF 
during the scaling matter and accelerated epochs. From 
Eqs. (14.3711 - (14.421) we obtain the perturbation equations 



X + 


V&q 

3xi 



(1 - C^)(lm 


(6.67) 


SN 3^ ~ 

C — -1" -Z^mVm , 

N 2 

x' = + i)x - C - ^, 

N 

S' = -3C'-IC^ (^"t) ’ 

= - (^H + V6qxi) Vm-^, 
m ^ 3»^L-2/c"(x + 0 

N 6 x? 


(6.59) 

(6.60) 
(6.61) 
(6.62) 
(6.63) 


where 1C = k/{dH). Taking the AA derivative of 
Eq. (j6.62|l and using other equations of motion, the ve¬ 
locity potential Vm obeys 


Let us derive solutions along which C and x stay nearly 
constant in the scaling matter era. Setting (C' ~ 0 and 
x' — 0 in Eqs. (16.5911 and (I6.60L respectively, we obtain 


. __ 3[1 - cg(l - y/ 6 qxi)]Clm ^ 

^ 2K?[2,clx\ +{I - cDen] 

t ^ 3[(1 + eg)(l - cl) + c^xi(3xi -I- '/&q)]Clm 

2iO[2,clxl + (1 — Ct)ei/] 

- _ clxi{?,xi - V&qen) s 

™ - t^[i^^xl +{I-cDen] 

5N _ _ 3c^ai(3xi - '/QqiH)^m j 

N 2 K?[?>clxl + {I - cDch] 


( 6 . 68 ) 

5m, (6.69) 

(6.70) 

(6.71) 


where, in the denominators of Eqs. (I6.68I) - (I6.71I) . we 
have neglected the terms without containing K?. The 
gravitational potentials and $, which are defined by 
Eq. (I4.43L obey 


i/"+Civ;; + 
v^g 


(1 C^)flrn ]*2 


2 clx\ 


JC^+C2 


Vm 


X- 






(6.64) 


d = ^5m , (6.72) 

2/C^ 

which is equivalent to Eq. (15.3011 with Eq. (15.281) . Since 
A 3 = 0 for gi = g 2 , the condition 2 HM‘^^/N jAg 
















































16 


used for the derivation of Eq. (15.301) is automatically sat¬ 
isfied. 

Substituting Eqs. (j6.68|) . (I6.71|) . and (j6.72l) into 
Eq. (j6.53L we obtain 


— 


2/C2 


-^ni,6n 


l + 2q^ 


(1 - c^){VGxi + 2q)(v^xi - 2qeH) 
2{3c?xf + (1 - cf)eH} 


(6.73) 


where xi = '/ 6 /[ 2 {fj, + q)] and in = 3^i/[2{fi + q)] for the 
fixed point (d) with w = 0. The effective gravitational 
coupling Geff during the scaling matter epoch reads 


Geff , ^ 2 I - Ct){V6xi+2q){^/6xl-2qiH) 

G ^ + 2{3c?i? + (1 - c2)e^} 

(6.74) 

where the term 2q^ arises in BD theories after the con¬ 
formal transformation to the EE [^. The last term on 
the rhs of Eq. (I6.74|) does not vanish for different from 
1. Since we are now considering the case where is 
constant, the variation of does not appear in the ex¬ 
pression of Eq. (16.741) . 

Eor the perturbations deep inside the Hubble radius, 
the rhs of Eq. ()5.19l) can be neglected relative to the Ihs 
of it. Then, during the scaling matter era, the matter 
perturbation obeys 


Provided that ^ q, there is a growing-mode solution 
to Eq. (16.751) . 




^^-‘ 2 q 
4(m + q) 


24(/r + g)^ . 

(/r-2g)2^ 


1 


(6.76) 

and g = (Geff/G)Hm- For q/g ^ 0 and g —>■ 1, the mat¬ 
ter density contrast evolves as dm oc d. In this case, the 
quantity {aH)^ 6 m, which appears in Eqs. (I6.68I) - (I6.71L 
remains constant, so the perturbations x, G Vm, and 
6 N/N do not vary in time. 

For q ^ 0 and Cj 1, the quantity g is different from 
1. As long as g <C 1 and g :s> q, the deviation of g 
from 1 is small, so the analytic formulas (I6.68|) - (I6.71I) 
are approximately valid in the scaling matter era. The 
perturbations x C start to vary after the Universe 
enters the epoch of cosmic acceleration, in which regime 
the analytic solutions (I6.68I) - (I6.71D are no longer valid. 

In Fig. m we plot the evolution of perturbations G Xj 
Vm, SN/N, and —(= l>) for the model parameters 
q = 0.1, c? = 0.5, Ai = 10, A 2 = 0.5, and U 2 /U 1 = 
10-6 (i.e., gi = 9.6, g 2 = 0.1, U 2 /U 1 = lO"®). The 
background initial conditions are chosen to start from the 
scaling matter fixed point (d), i.e., xi = •\/6/[2(Ati -|- g)] 
and X 2 = \/(3 -I- 2q'^ + 2qgi)/[2{gi + g)^]. We choose 



a 


Figure 1. Evolution of the perturbations G Xi Im, 5N/N, 
and —'I' in the EF versus the scale factor a for the model 
parameters q — 0.1, c? = 0.5, Ai = 10, A 2 = 0.5, and V 2 /V 1 = 
10“®. We choose the initial values of xi and £2 to coincide 
with those of the scaling fixed point (d) with w = 0 and gi = 
q 2 = 0.1. The initial values of perturbations are chosen to be 
^ = 6.8407 X 10"®, X = -2.7513 x 10“®, Vm = 3.5537 x 10“®, 
d = 2.5618 X 10"®, and /C = 30 at a = 2.0488 x 10“®, under 
which the special solution (16.671) is the dominant contribution 
to Vm- 


the initial value of Vm close to the special solution (16.671) 
with ~ 0 and x' ~ 0 for the normalized wave number 
1C = 30. 

All the perturbations shown in Fig. [1] stay nearly con¬ 
stant during the scaling matter era. We also confirmed 
that the analytic formulas (I6.68l) - (|6.72l) are in good agree¬ 
ment with numerically integrated solutions in the scal¬ 
ing regime. For sub-horizon perturbations (/C 1), 

the choice of different wave numbers k only modifies 
the amplitudes of perturbations Xi G Vm \ and SN/N. 
There is a simple relation SN/N = — {sClm/2)Vm'^ from 
Eqs. (16.701) and ()6.71l) . The perturbations SN/N and 
Vm are suppressed relative to x and ( because of the 
conditions xi 1 and g <C 1. The Universe finally en¬ 
ters the epoch of cosmic acceleration driven by the fixed 
point (c) with ^2 = A 2 — 4g. As we see in Fig. [U the 
perturbations start to vary after the onset of cosmic ac¬ 
celeration. 

Numerically, we also compute the gravitational poten¬ 
tial and the resulting effective gravitational coupling 
Geff from the definition (16.521) . In Fig. [5] the evolution of 
Geff/G is plotted for four different values of g and c^. We 
confirmed that the analytic estimation (16.741) of Geff is in 
good agreement with the numerical result in the scaling 
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Figure 2. Evolution of the effective gravitational coupling Gefr 
normalized by the gravitational constant G in the EF for the 
four cases: (i) q = 0, Ct = 0.5, (ii) q = 0.05, c? = 0.5, (iii) 
q = 0.1, cj = 0.5, and (iv) q = 0.1, cf = 0.99, with the model 
parameters Ai = 10, A 2 = 0.5, and V 2 /V 1 = 10“®. The initial 
conditions are chosen in the same way as those explained in 
the caption of Fig. [T] 


matter regime. The growth of GeS starts to occur after 
the scaling matter era. 

If we compare the cases (ii) and (iii) with the case (i) 
in Fig. [21 we find that the existence of coupling q leads to 
the value of Ges/G smaller than that for q = 0. Provided 
Cj is not close to 1, the analytic formula (I6.74|) implies 
that Ges/G approaches 1 in the limit where q^ xi and 
xi <C 1. Hence, for a given value of cl different from 1, 
Geff/G tends to decrease with increasing q. For larger 
g, the variation of GeS/G occurs at a later cosmological 
epoch. 

In the limit that 1, Eq. ()6.74|) reduces to the value 

I 


Gefi/G —>■ 1 + 2g^. The case (iv) in Fig.|2]is close to such 
an example, in which case the variation of Geff/G is small 
even after the onset of cosmic acceleration. This property 
can be clearly distinguished from the model with g = 0 
and cl ^ I. 

If we choose initial conditions where Vm is not close to 
the special solution the homogenous solution Vm'^ 
gives rise to the oscillation of Vm- This oscillation con¬ 
tinues by the time when the amplitude of V,^'^ decreases 

^ fs) 

sufficiently relative to that of Vm ■ This situation is anal¬ 
ogous to what was found for the case g = 0 [s^ . 

2. Jordan frame 


The evolution of perturbations in the EF can translate 
to that in the JF by using the correspondence 


X = 




SN = 


1 -I- OJh' 
SN 


N 


c = c, 


Vm 


Vm 

1 -I- UJh ’ 


S = 5, ^ = ^ 

Ct 


where/C = k/{aH) and w// = — •\/6gii/(l + '\/6gii). The 
gauge-invariant matter perturbation Sm is related to Sm, 
as 


Sm = L + -^^Vm . (6.78) 

l + UJH 

At the background level, we also have the relations ()6.38l) 
and (16.401) . Using the analytic solutions (I6.68I) - (I6.71D in 
the EF during the scaling matter era, the perturbations 
X, G Vm, and SN in the JF can be expressed in terms of 
/C, rim, cl, xi, and ch- 

The gravitational potentials dt and $ are known from 
4>, and x by using Eqs. (16.571) and (16.581) . Substituting 
the solutions (I6.68P and (16.721) into Eqs. (I6.57|) and (16.581) 
in the scaling matter regime and employing the approx¬ 
imation Sm — Sm for the perturbations deep inside the 
Hubble radius, it follows that 


‘S[{l-cl)eH{l + ‘ 2 .ujH) + ‘ix\+ujlj\ ^ 

^ - 2/C2[(I - cl){eH + cc^)(I + u:h) + icW^ ’ 

dj ^ 3 [c^(l -|- 2 u}h) + (I + Wi/){I + tH 2 .UIH — Cj (2 + tH Swi/)} -|- Sclxf] 

2 c^/C 2 [(l - cI)(€h + UJh )(1 + OJh) + 3 clxj] ™ ’ 

[(1 - cl){l + eg) -f 2 ujh ][1 -cl+ ojHjl - 2c^)] 
cl[{l - cDenil+ 2 ojh) + 3 x 1 + uj%] 

^ {1 — cI)€h{ 1 ^ 2 ujh) 3 x 1 + ojf^ G 

“ (1 - cl){eH + ojh){1 + OJh) + 3clxj Fi ’ 


(6.79) 

(6.80) 
(6.81) 
(6.82) 


where ojh 


—'/Gqxi and xi = •\/6/(2Ai) ^ \'Xh\ for g <C I. The perturbation x induces the anisotropic stress 
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Figure 3. The gravitational potentials —Vlt, and —3>eff = 
($ — ^)/2 in the JF versus 1 + z = Ija for the same model 
parameters and initial conditions as those given in Fig.[T] The 
present epoch (z = 0) is identified by the condition flrn = 0.3. 
The vertical dot-dashed line represents the onset at which the 
cosmic acceleration (a > 0) sets in (z ~ 0.63). Numerically, 
we integrate the background and perturbation equations in 
the EF and then compute 'F and <1> by using the transforma¬ 
tion laws (16.571) and (16.581) . 


in the JF, so the anisotropy parameter 77 is different from 

1 . 


In the limit —)► 1, the formulas (I6.79I) - (I6.82I) give ~ 
-3(1 + 2g^)Q^S,n/(2/C^), $ 3(1 - 2q^)nm5ml{210). 

and 


77 ~ 


1 - 2g2 
l + 2g2 


Geff — (1 + 2 q^)— . 


(6.83) 

(6.84) 


In the limit g —>■ 0 (i.e., ojh —>■ 0), we use the approxi¬ 
mations ~ 3/2 and flm — 1 during the scaling regime 
and eliminate the term xl on account of Eq. (EH). This 
process leads to 


77 ~ 1 -h 


Geff ^ 1 


5(l-cg)(c2-c,^) 
3c?(l + c2 - cf) 


1-c? 


G, 


(6.85) 

( 6 . 86 ) 


which match those derived in Ref. without referring 
to the EF. Since Cg can be much greater than 1 for 1, 
the parameter 77 exhibits the large deviation from 1. In 
this case, Geff is slightly larger than G. 

If yf 1 and g yf 0, the difference between the gravita¬ 
tional potentials — and depends on the magnitudes of 
the terms 1 — and Wij. Provided |1 — Cj ] {\ujh\.x\}, 


Figure 4. Evolution of f{z)as{z) versus the redshift 2 in the 
JF for Ai = 10, A 2 = 0.5, and V 2 /V 1 — 10“®. From the top 
to the bottom, each curve corresponds to (i) (7 = 0, c? = 0.5, 
(ii) q = 0.05, Ct — 0.5, (iii) q = 0.1, cj = 0.5, and (iv) 
q = 0.1, Ct = 0.99, respectively. The initial conditions are 
chosen in the similar way to those explained in the caption of 
Fig.ID For a given value of c? (yf 1), the growth rate of mat¬ 
ter perturbations tends to be smaller with increasing q due to 
the decrease of Geff. The black points with error bars corre¬ 
spond to the data from the recent observations of f(z)crs(z) 
by 2dFGRS [i^ 6 dFGRS [lH, WWleZ [H, SDSSLRG 
BOSSGMASS0, and VIPERS [HI surveys. 


the parameter 77 of Eq. (16.811) is close to the value (16.831) . 
whereas, for |1 — c^\ ^ {\ojh\,xI}, rj is close to the value 

EISD- 

In Fig. [3] we illustrate the evolution of —T, $, and 
—$eff = (4> — lE')/2 versus the redshift z in the JF for 
Cj = 0.5, q = 0.1, and Ai = 10. During the scaling matter 
epoch, the gravitational potentials stay nearly constant. 
Since |1 —| {\ijjH\.x\} for the model parameters used 

in the simulation of Fig. O we have that 77 ~ 2.6 from 
Eq. (16.851) . This is in good agreement with the numer¬ 
ical result of Fig. [3] in the scaling matter regime. The 
gravitational potentials start to decrease after the onset 
of cosmic acceleration, but —df shows temporal growth 
in the future because of the increase of GeS (see the case 
(iii) in Fig. El). Finally, the Universe enters the attractor 
regime in which all of —4', 4), and —$eff decrease in a 
similar way. 

In Fig. m we plot the evolution of f{z)as{z) in the 
low-redshift regime for four different values of q and 
(corresponding to those used in Fig. [3]), where f{z) = 
Sm/{HSm) and as{z) is the rms amplitude of over-density 
at the comoving 8 h~^ Mpc scale {h is the normalized 
Hubble parameter Hq = 100 h kmsec“^ Mpc“^). This 
quantity is associated with the observations of red-space 
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distortions in galaxy clusterings. Note that the growth 
rate f{z)a^{z) can be also measured by using only the 
peculiar motions of galaxies in the low redshift [85|, [Sq . 
The initial condition of 6 m is chosen such that its am¬ 
plitude today is equivalent to fT8(0) = 0.82 [s^. During 
the scaling matter era, we have numerically checked that 
Geff is well described by Eq. (j6.82l) . When q = 0 and 
Cj = 0.5 we have Ges — 1.03G from Eq. (16.861) . whereas, 
for q = 0.1 and ~ 1, GeS — 1.02G/Fi from Eq. (16.841) . 

For larger q with a given value of 1), the onset 

of growth of Geff occurs at later cosmological epochs. In 
fact, Fig.lHshows that the values of f{z)as{z) for = 0.5 
tend to be smaller with increasing q in the low-redshift 
regime. For Cj close to 1, the variation of Geff is small 
by the present epoch (see Fig. [2]). In this case, Geff is 
approximately given by Eq. (|6.84l) even around today. 

In Fig. |4] we also plot the recent data of f{z)as{z) 
constrained by the redshift-space distortion (RSD) mea¬ 
surements. Using the bound on (T8(0) from the Planck 
data [13 , the RSD data tend to favor the growth rate of 
6 m lower than that predicted by GR [s^. In Fig. |4l such 
a property can be also observed in our model where Geff 
is slightly larger than G. The recent 6dF galaxy surveys 
using only the peculiar motions of galaxies provided the 
constraint /(0)tT8(0) = 0.415 ± 0.065, which is consistent 
with the four cases shown in Fig. 01 It remains to see 
how future RSD and peculiar velocity measurements will 
pin down the error bars of f(z)as{z). 

We have thus clarified the evolution of observables as¬ 
sociated with the measurements of CMB, redshift-space 
distortions, and weak lensing by transforming back from 
the EF to the JF. Since the EFTCAMB code [s^ for 
modified gravity theories is written in the JF, our re¬ 
sults in this section can be used to place observational 
constraints on the model (15.221) with the functions (16.11) . 
Since the staring point of coupled scalar-field models (in¬ 
cluding coupled quintessence chameleons @ , and 
disformally coupled models |31|-|34|1 is usually assumed 
to be the EF, our analysis in the EF is also useful to con¬ 
strain coupled dark energy models with Cj different from 
1. We leave observational constraints on such models for 
a future work. 


VII. CONCLUSIONS 

In the presence of matter, we have studied cosmological 
disformal transformations in a generalized class of Horn- 
deski theories (GLPV theories). In these theories there 
is one propagating scalar degree of freedom (p coupled to 
the metric g^i, in the JF on the flat FLRW background. 
Even if matter is minimally coupled to gravity in the 
JF, the matter sector feels the modification of gravity 
through the change of gravitational potentials mediated 
by the field </>. 

The structure of the Lagrangian in GLPV theories, 
which is given by Eq. (1^ in the unitary gauge, is pre¬ 
served under the disformal transformation (HH), while 


the matter Lagrangian contains a coupling with the field 
(j) and its derivatives in the transformed frame. Thus, 
the matter-scalar interaction becomes explicit after the 
disformal transformation. In Sec. m we clarified how 
the energy-momentum tensor of matter and associated 
background/perturbed quantities are mapped under the 
disformal transformation. 

In Sec. m we have derived the background and lin¬ 
ear perturbation equations of motion in both the JF and 
the transformed frame. In the transformed frame, the 
coupling Q in Eq. (14.361) arises for the matter continuity 
Eq. (14.351) at the background level. The matter perturba¬ 
tion Eq. ()4.13l) in the JF is also transformed to the more 
involved Eq. (14.411) due to the matter-scalar interaction, 
while the structure of other perturbation equations is not 
subject to change. 

In Sec.|V]we discussed the transformation from the JF 
to the EF in which the second-order action of tensor per¬ 
turbations is of the same form as that in GR. Under the 
choice (15.51) of the factors H and P, the Bardeen poten¬ 
tials 'L and $ in the EF obey the “de-mixed” relation 
(15.161) . If the action in the EF belongs to a class of Horn- 
deski theories (dn = 0), there is no anisotropic stress 
between T and $. 

The non-relativistic matter density contrast 6 m obeys 
the differential Eq. (15.191) in the EF, where 'Lg is the 
effective gravitational potential given by Eq. (I5(20p . In 
the EF, it becomes transparent that the variations of 
H and Ct as well as the deviation of Cf from 1 lead to 
the modification of gravitational interactions with matter 
perturbations. The gravitational potential T in the JF 
is simply related to Tg, as T = C^'Lg. 

In Sec. EH we proposed a concrete model of dark en¬ 
ergy in which the coupling between matter and the scalar 
degree of freedom tp is manifest after the disformal trans¬ 
formation to the EF. For the field potential (16. 8p with 
Ai > 10 and A 2 ;< 1, there exist scaling solutions corre¬ 
sponding to radiation and matter eras followed by an at¬ 
tractor with the cosmic acceleration. At the background 
level, the disformal transformation to the EF gives rise 
to the term q 2 associated with the function but we 
showed that the stability of fixed points is independent of 
< 72 . This reflects the fact that the background equations 
in the JF do not contain the function B^. 

We also studied the evolution of linear cosmological 
perturbations from the matter era to today for the case 
Qi = 92 - In the EF we derived the second-order equation 
of the velocity potential Vm and identified the special so- 
lution Vm on scales deep inside the Hubble radius. For 
the initial conditions satisfying \Vm\, we ob¬ 

tained analytic solutions of perturbations where C and x 
stay nearly constant. On using these solutions, we de¬ 
rived the effective gravitational potential ^g of the form 
(16.731) during the scaling matter era. The coupling q and 
the deviation of from 1 lead to the gravitational cou¬ 
pling Geff modified from that in GR. 

Once the evolution of perturbations is known in the 
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EF, it is straightforward to transform it back to that 
in the JF by using the correspondence (16.571) - (16.581) and 
(I6.77I) - (I6.78I) . While $ = —'k in the EF for the La- 
grangian (15.271) . the field x generates the anisotropic 
stress in the JF such that r] = —$/'k ^ 1. For sub¬ 
horizon perturbations the effective gravitational coupling 
Geff in the JF is related to Geff in the EF as Ges — 
Ges/Fi, so the growth rate of matter perturbations in 
the JF is known accordingly. 


We have analytically estimated rj and Geff during the 
scaling matter era and confirmed that they are in good 
agreement with numerical results before the onset of cos¬ 
mic acceleration. We also numerically computed the evo¬ 
lution of f(z)as{z) by transforming back from the EF to 
the JF. As we see in Fig. HI it is possible to distinguish 
between the models with different values of q and ob- 
servationally. 


We have thus shown that the disformal transformation 
is useful for understanding gravitational interactions with 
matter mediated by the scalar field. After transforming 
to the EF, the background and perturbation dynamics 
in the JF are readily known by using the correspondence 
of physical quantities between the two frames. We can 
apply our results to observational constraints on dark en¬ 
ergy models in the framework of GLPV theories. More¬ 
over, our analysis in the EF is useful for constraining 
coupled dark energy models in which the starting point 
is the EF Lagrangian with matter-scalar couplings. 


While we focused on the cosmological set up, it is of 
interest to extend the disformal transformation to gen¬ 
eral space-time including the spherically symmetric back¬ 
ground. This should help us to understand the nature of 
matter-scalar couplings in local regions of the Universe. 
In particular, the derivation of the effective gravitational 
coupling around a compact body (like the Sun) will be 
important to place constraints on theories beyond Horn- 
deski from local gravity experiments. 
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Appendix A: Correspondence between the two 
frames 


We show relations for the quantities connected through 
the disformal transformation. The background quantities 
are transformed as 
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For the quantities associated with perturbations, we have 
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